Nonlinear shear-current dynamo and magnetic helicity transport in sheared 

turbulence 
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The nonlinear mean-field dynamo due to a shear-current effect in a nonhelical homogeneous 
turbulence with a mean velocity shear is discussed. The transport of magnetic helicity as a dynamical 
nonlinearity is taken into account. The shear-current effect is associated with the WxJ term in 
the mean electromotive force, where W is the mean vorticity due to the large-scale shear motions 
and J is the mean electric current. This effect causes the generation of large-scale magnetic field 
in a turbulence with large hydrodynamic and magnetic Reynolds numbers. The dynamo action 
due to the shear-current effect depends on the spatial scaling of the correlation time r(fe) of the 
background turbulence, where k is the wave number. For Kolmogorov scaling, r(A;) oc k~^^^, the 
dynamo instability occurs, while when r(fc) cx (small hydrodynamic and magnetic Reynolds 
numbers) there is no the dynamo action in a sheared nonhelical turbulence. The magnetic helicity 
flux strongly affects the magnetic field dynamics in the nonlinear stage of the dynamo action. 
Numerical solutions of the nonlinear mean-field dynamo equations which take into account the 
shear-current effect, show that if the magnetic helicity flux is not small, the saturated level of the 
mean magnetic field is of the order of the equipartition field determined by the turbulent kinetic 
energy. Turbulence with a large-scale velocity shear is a universal feature in astrophysics, and 
the obtained results can be important for elucidation of origin of the large-scale magnetic fields in 
astrophysical sheared turbulence. 
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1. Introduction 

It is generally believed that one of the main reasons for 
generation of large-scale magnetic fields in turbulent flow 
is the a effect (see, e.g., books and reviews by MofFatt 
1978, Parker 1979, Krause and Radler 1980, Zeldovich 
et al. 1983, Ruzmaikin et al. 1988, Stix 1989, Roberts 
and Soward 1992, Ossendrijver 2003, Brandenburg and 
Subramanian 2005a). However, the a effect caused by 
the helical random motions of conducting fluid, requires 
rotating inhomogeneous or density stratifled turbulence. 

In a turbulence with a large-scale velocity shear and 
high hydrodynamic and magnetic Reynolds numbers 
there is a possibility for a mean-field dynamo (see Ro- 
gachevskii and Kleeorin 2003, 2004). Turbulence with a 
large-scale velocity shear is a universal feature in astro- 
physical plasmas. The large-scale velocity shear creates 
anisotropic turbulence with a nonzero background mean 
vorticity W. This can cause the WxJ effect (or the 
shear-current effect), which creates the mean electric cur- 
rent along the original mean magnetic field and produces 
the large-scale dynamo even in a nonrotating and non- 
helical homogeneous turbulence (see Rogachevskii and 
Kleeorin 2003, 2004). Here J is the mean electric cur- 
rent. 

The mean-field dynamo instability is saturated by the 
nonlinear effects. A dynamical nonlinearity in the mean- 
field dynamo which determines the evolution of small- 
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scale magnetic helicity, is of a great importance due to 
the conservation law for the total (large and small scales) 
magnetic helicity in turbulence with very large magnetic 
Reynolds numbers (see, e.g., Kleeorin and Ruzmaikin 
1982, Gruzinov and Diamond 1994, 1996, Kleeorin et 
al. 1995, 2000, 2002, 2003a, 2003b, Kleeorin and Ro- 
gachevskii 1999, Blackman and Field 2000, Vishniac and 
Cho 2001, Blackman and Brandenburg 2002, Branden- 
burg and Subramanian 2005a, Zhang et al. 2006). On 
the other hand, the effect of the mean magnetic field on 
the motion of fluid and on the cross-helicity can cause 
quenching of the mean electromotive force which deter- 
mines an algebraic nonlinearity. The combined effect of 
the dynamic and algebraic nonlinearities saturates the 
growth of the mean magnetic fleld. 

The mean-field dynamo is essentially nonlinear due to 
the evolution of the small-scale magnetic helicity (Gruzi- 
nov and Diamond 1994, 1996). In particular, even for 
very small mean magnetic field the magnetic a effect is 
not small. This is a reason why we have to take into 
account the dynamical nonlinearity in the mean-field dy- 
namo. When we ignore the dynamical nonlinearity due 
to evolution of small-scale magnetic helicity and take 
into account only algebraic nonlinearity caused by the 
nonlinear shear-current effect, we obtain the saturated 
level of mean magnetic field which is in several times 
larger than the equipartition field determined by the tur- 
bulent kinetic energy (Rogachevskii et al. 2006). This 
result can be important in view of astrophysical appli- 
cations whereby the super-equipartition large-scale mag- 
netic fields are observed, e.g., in the outer parts of a few 
galaxies (see Beck 2004, 2005). Note that it is a prob- 
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FIG. 1: Mechanism for the a effect: (a). Interaction between 
uniform original mean magnetic field and inhomogencous ro- 
tating turbulence; (b). The deformations of the original mag- 
netic field lines are caused by the upward and downward tur- 
bulent eddies; (c). Formation of the mean electric current 
£" = aB oc —Iq (f2-A")B oc J opposite to the original mean 
magnetic field (for $7- A" > 0). Here Ju and Jd are the electric 
currents caused by the deformations of the original magnetic 
field lines by the upward and downward turbulent eddies, re- 
spectively, and |Ju| > |Jd|. 



FIG. 2: Mechanism for the WxJ effect: (a). Interaction 
between nonuniform original mean magnetic field and homo- 
geneous sheared turbulence; (b). The deformations of the 
original magnetic field lines are caused by the upward and 

downward turbulent eddies; (c) . Formation of the mean elec- 
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the original mean magnetic field (for W-A^ > 0). Here 
and Jd are the electric currents caused by the deformations of 
the original magnetic field lines by the upward and downward 
turbulent eddies, respectively, and |J„| > |Jd|. 



lorn to roach a supcr-cquipartition level of the large-scale 
magnetic field in the afl dynamo. 

The goal of this study is to investigate the nonlinear 
mean-field dynamo due to the shear current effect. In 
this study we have taken into account the dynamic non- 
linearity caused by the evolution of the small-scale mag- 
netic helicity. This paper is organized as follows. In 
section 2 wc elucidate the physics of the shear-current ef- 
fect. In section 3 we consider kinematic dynamo problem 
due to this effect. In section 4 we describe the results of 
numerical solutions of the nonlinear mean-field dynamo 
equations which take into account the shear-current ef- 
fect, the algebraic and dynamic nonlinearities. Finally, 
the discussion and conclusions are given in section 5. In 
Appendixes A we compare the problem of generation of 
the mean magnetic field in a turbulence with a large- 
scale velocity shear with that of the generation of mean 
vorticity in a sheared turbulence. 



2. The physics of the shear-current effect 

In this section we elucidate the mechanism of genera- 
tion of large-scale magnetic field due to the shear-current 
effect. To this end wc first discuss the physics of the a 
effect (see, e.g., Moffatt 1978, Parker 1979, Krause and 
Radler 1980, Zeldovich et al. 1983, Ruzmaikin et al. 
1988). The a term in the mean electromotive force E = 
(u X b) in a rotating inhomogencous turbulence can be 
written in the form = aB oc —Iq (r2-A")B (see, e.g., 
Krause and Radler 1980, Radler et al. 2003), where u 
and b are fluctuations of the velocity and magnetic field, 
respectively, angular brackets denote ensemble averaging, 
f2 is the angular velocity, the vector A" = V(u^)/(u^) 
determines the inhomogeneity of the turbulence, B is the 
mean magnetic field and /q is the maximum scale of tur- 
bulent motions (the integral turbulent scale). The a ef- 
fect is caused by the kinetic helicity Xu oc r}^ (O-A") in 
an inhomogencous rotating turbulence, where -q^ cx Iq uq 
is the turbulent magnetic diffusion and uq is the charac- 



FIG. 3: Mechanism for shear-induced generation of pertur- 
bations of the mean magnetic field By by sheared stretching 
of the field Bx- This effect is determined by the first term 
(oc SBx) in RHS of Eq. (|5J, and it is similar to the differen- 
tial rotation because Vx(UxB) = SB^Sy. 



teristic turbulent velocity in the maximum scale of turbu- 
lent motions Iq. The deformations of the magnetic field 
lines are caused by upward and downward rotating tur- 
bulent eddies (see figure 1). The inhomogeneity of the 
turbulence breaks a symmetry between the upward and 
downward eddies. Therefore, the total effect of these ed- 
dies on the mean magnetic field does not vanish, and it 
creates the mean electric current along the original mean 
magnetic field due to the a effect. 

The large-scale magnetic field can be generated even 
in a nonrotating and nonhelical turbulence with a mean 
velocity shear due to the shear-current effect (see Ro- 
gachevskii and Kleeorin 2003, 2004). This effect is related 
to the WxJ term in the mean electromotive force, and 
it can be written in the form 5* oc — ZqWx(VxB) oc 

(W-A^)B, where the mean vorticity W — VxU is 
caused by the mean velocity shear and = VB^/2B^ 
determines the inhomogeneity of the mean original mag- 
netic field. In a sheared turbulence the inhomogeneity of 
the original mean magnetic field breaks a symmetry be- 
tween the influence of upward and downward turbulent 
eddies on the mean magnetic field. The deformations of 
the original magnetic field lines in the WxJ effect are 
caused by the upward and downward turbulent eddies. 
This creates the mean electric current along the mean 
magnetic field and produces the magnetic dynamo (see 
figure 2). 

In order to demonstrate how the shear-current dynamo 
operates, let us consider a homogeneous turbulence with 
a mean velocity shear, U = (0, Sx, 0) and W = (0, 0, S). 
Let us assume that the mean magnetic field has a simple 
form B = {Bx{z), By{z),0). The mean magnetic field in 
the kinematic approximation is determined by 
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dt 



FIG. 4: Mechanism for shear-current generation of perturba- 
tions of the mean magnetic field Bx from the inhomogeneous 
magnetic field By . This effect is determined by the first term 
in RHS of Eq. Q. 
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(Rogachevskii and Kleeorin 2003, 2004), where B'^ = 
d^Bi/dz^, T]^ is the coefficient of turbulent magnetic dif- 
fusion and the dimensionless parameter determines 
the WxJ effect (see Eq. (O below). The first term 
oc SB, in the right hand side of Eq. ||2Jl determines 
the stretching of the magnetic field Bx by the shear mo- 
tions, which produces the field By (see figure 3). On the 
other hand, the interaction of the non-uniform magnetic 
field By with the background vorticity W (caused by the 
large-scale shear) produces the electric current along the 
field By. This implies generation the field component Bx 
(see figure 4) due to the WxJ effect, which is determined 
by the first term in the right hand side of Eq. . This 
causes the dynamo instability. 



3. Kinematic dynamo due to the shear-current effect 

Let us consider the kinematic dynamo due to the shear- 
current effect. In this study we have derived a more gen- 
eral form of the parameter for arbitrary scaling of 
the correlation time T(fc) of turbulent velocity field. The 
generalized form of the parameter defining the shear- 
current effect, is derived using Eq. (A44) given in the 
paper by Rogachevskii and Kleeorin (2004) . The param- 
eter (Tg entering in Eq. |^ is given by 
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where 



T'(fc) k T{k) E{k) dk , T^{k) E{k) dk , 



-a^Sll B'' + 7j, B'^ 



(1) 



(4) 

E{k) is the turbulent kinetic energy spectrum, T(fc) is 
the scale-dependent correlation time of turbulent velocity 
field, the parameter e = Em Im/Ey Iq, Em and Ey are the 
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magnetic and kinetic energies per unit mass in the back- 
ground turbulence (with a zero mean magnetic field), Im 
is the scale of localization of the magnetic fluctuations 
generated by the small-scale dynamo in the background 
turbulence and T'{k) = dr/dk. Equations lO and Q are 
written in dimensionless form, where the turbulent time 
T(fc) is measured in the units of tq = Io/uq, the wave 
number k is measured in the units of Iq^, and the tur- 
bulent kinetic energy spectrum E{k) is measured in the 
units of Uq Iq. 

The solution of equations Q and ((SJ we seek for in 
the form oc exp(7t -|- iKzz), where the growth rate, 7, 
of the mean magnetic field due to the magnetic dynamo 
instability is given by 

J = Slo^/a^ K;, ~ 1]^ . (5) 

The necessary condition for the magnetic dynamo insta- 
bility is > 0. The dynamo action due to the shear- 
current effect depends strongly on the spatial scaling of 
the correlation time T{k) of the turbulent velocity field. 
In particular, when T(fc) oc fc^^, the ratio I/Iq — — /i, 
and the criterion for the dynamo instability reads 

1 - /i + 3 e > . (6) 

For example, when T(fc) cx fc~^/^ (the Kolmogorov scal- 
ing), the parameter cr^ = (4/135) (1 4-9e). This case was 
considered by Rogachevskii and Kleeorin (2003; 2004). 
When e = (there are no magnetic fluctuations in the 
background turbulence due to the small-scale dynamo), 
the shear-current dynamo occurs for ^ < 1. The bound- 
ary 11 — 1 corresponds to the spatial scaling of the 
correlation time T(fc) cx fc~^. For the Kolmogorov's 
type turbulence (i.e., for a turbulence with a constant 
energy flux over the spectrum), the energy spectrum 
E{k) — —dr/dk oc k~^. This implies that the veloc- 
ity is dominated by the large scales more strongly than 
for the turbulence with a purely Kolmogorov spectrum 
E{k) oc fc-5/3. 

For small hydrodynamic and magnetic Reynolds num- 
bers, the turbulent time T{k) oc l/(i^fc^) or T(fc) oc 
1/(77/0^) depending on the magnetic Prandtl number, i.e., 
T(fc) oc Thencj^ = (4/o/15) (-1-1-3 e), where ^ 2, 
v is the kinematic viscosity and rj is the magnetic diffu- 
sivity due to electrical conductivity of the fluid. When 
e = the parameter < 0, and there is no dynamo 
action due to the shear-current effect in agreement with 
the recent studies by Radler and Stepanov (2006) and 
Riidiger and Kichatinov (2006). They have not found 
the dynamo action in nonrotating and nonhelical shear 
flows with e = in the framework of the second order 
correlation approximation (SOCA). This approximation 
is valid for small hydrodynamic Reynolds numbers. Even 
in a highly conductivity limit (large magnetic Reynolds 
numbers), SOCA can be valid only for small Strouhal 
numbers, while for large hydrodynamic Reynolds num- 
bers (i.e., for fully developed turbulence), the Strouhal 
number is 1. When e > 1/3 the mean magnetic field can 



be generated due to the shear-current effect even for small 
hydrodynamic and magnetic Reynolds numbers. How- 
ever, the latter case seems to be not realistic. 

The effect of shear on the mean electromotive force and 
shear-current effect have been studied by Rogachevskii 
and Kleeorin (2003, 2004) for large hydrodynamic and 
magnetic Reynolds numbers using two different ap- 
proaches: the spectral r approximation (the third-order 
closure procedure) and the stochastic calculus, i.e., the 
Feynman-Kac path integral representation of the solution 
of the induction equation and Cameron-Martin-Girsanov 
theorem. 

Note that Ruderman and Ruzmaikin (1984) formally 
constructed an example of an exponentially growing mag- 
netic field in a fluid with shear and a homogeneous 
anisotropic magnetic diffusivity. An essential condition 
for generation is that the vector defining the anisotropy 
in this phenomenological model must be non-parallel 
and non-perpendicular to the velocity. However, equa- 
tions and Q arc different from those given by Ruder- 
man and Ruzmaikin (1984) because they have not study 
the effect of shear on the mean electromotive force. The 
first (although incorrect) attempt to determine the effect 
of shear on the mean electromotive force has been made 
by Urpin (1999a, 1999b) in the framework of SOCA. 

In order to study the kinematic dynamo due to the 
shear current effect, let us rewrite equations ^ and Q 
for the mean magnetic field in the dimensionless form 

where the mean magnetic field is B = By(t, z) By + 
S~'^Vx[A{z)ey], i.e., B^{t,z) = -S'^M'(z), the param- 
eter S"* = S L'^/r/j, is the dimensionless shear number and 
D — {Iq S^,/ L)^ a g is the dynamo number. We consider 
the following boundary conditions for a layer of the thick- 
ness 2L in the z direction: 

By{t,\z\ = l)=Q- ^'(t,|z| = l)-0, (9) 

i.e., B(i, \z\ = 1) = 0. In dimensionless equations 
and ||SJ| the length is measured in units of L, the time 
is measured in units of the turbulent magnetic diffusion 
time L'^/rjj,, the mean magnetic field B is measured in 
units of the equipartition field Bcq — ^/AiFpuo deter- 
mined by the turbulent kinetic energy and the turbu- 
lent magnetic diffusion coefficient is measured in units of 
the characteristic value of t]^ = IqUq/3. The solution of 
equations {T)) and |(SJ) reads 

By{t,z) ^ Bo e-Kp{'yt) cos{KzZ + if) , (10) 

Bx{t,z) = ^ ^/ctJBo exp(7t) cos(X2Z + (^)(11) 

The growth rate of the mean magnetic field in the dimen- 
sionless form is given by 7 = \/D — . The wave 
vector Kz is measured in units of i^-'^ and the growth 
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rate 7 is measured in units of the inverse turbulent mag- 
netic diffusion time rj^/L^. 

For the symmetric mode the angle (p = irn, the wave 
number Kz = (7r/2)(2m + 1) and the mean magnetic 
field is generated when the dynamo number D > I?cr — 
(7r^/4)(2rn + 1)^, where n,m — 0, f,2,.... For this 
mode the mean magnetic field is symmetric relative to 
the middle plane z = 0. For the antisymmetric mode 
the angle ip = (7r/2) (2n + 1) with n = 0,1,2,..., the 
wave number Kz = tt m and the magnetic field is gen- 
erated when the dynamo number D > D^r — tt^ , 
where m = 1,2,3,.... Note that for the shear-current 
dynamo, the ratio of the field components B^/By is 
small [i.e., B^/By ^ jlg/L) Kz ^fa^ <C 1 when Kz is 
not very large], see equations (|10|l and Hll|) . This fea- 
ture is similar to that for the a£l dynamo, whereby the 
poloidal component of the mean magnetic field is much 
smaller than the toroidal field. The maximum growth 
rate of the mean magnetic field, 7inax = £'/4, is at- 
tained at Kz = Kra — This corresponds to the 
characteristic scale of the mean magnetic field variations 
LB = 2TiL/K^ = {ATT/y/D)L. 



4. Nonlinear dynamo due to the shear-current effect 

Kinematic dynamo models predict a field that grows 
without limit, and they give no estimate of the magni- 
tude for the generated magnetic field. In order to find the 
magnitude of the field, the nonlinear effects which limit 
the field growth must be taken into account. The nonlin- 
ear theory of the shear-current effect was developed by 
Rogachevskii and Kleeorin (2004, 2006). 



4.1. The algebraic nonlinearity 

First, let us start with the algebraic nonlinearity which 
is determined by the effects of the mean magnetic field 
on the motion of fluid and on the cross-helicity. These 
effects cause quenching of the mean electromotive force. 

Below we outline the procedure of the derivation of 
the equations for the nonlinear coefficients defining the 
mean electromotive force in a homogeneous turbulence 
with a mean velocity shear (see for details, Rogachevskii 
and Kleeorin 2004). We use the momentum equation and 
the induction equation for the turbulent fields written in 
a Fourier space. We derive equations for the correla- 
tion functions of the velocity field fij — {uiUj), the mag- 
netic field hij — (bibj) and the cross-helicity gij = (btUj), 
where the angular brackets denote ensemble averaging. 
We split the tensors fij, hij and gij into nonhelical, hij, 

and helical, hl^ , parts. The helical part of the tensor 

hij for magnetic fluctuations depends on the magnetic 
helicity, and it is determined by the dynamic equation 
which follows from the magnetic helicity conservation ar- 
guments (see Section 4.2). Then we split the nonheli- 



cal parts of the correlation functions fij , hij and gij into 
symmetric and antisymmetric tensors with respect to the 
wave vector k. 

The second-moment equations include the first-order 
spatial differential operators JV applied to the third-order 
moments M'^^^^\ A problem arises how to close the sys- 
tem, i.e., how to express the set of the third-order terms 
MM'^^^^^ through the lower moments M^-^-^^ (see, e.g., 
Orszag 1970; Monin and Yaglom 1975; McComb 1990). 
Various approximate methods have been proposed in or- 
der to solve it. A widely used spectral r approximation 
(see, e.g., Orszag 1970, Pouquet et al. 1976, Kleeorin 
et al. 1990, Kleeorin et al. 1996, Blackman and Field 
2002, Rogachevskii and Kleeorin 2004, Brandenburg et 
al. 2004, Brandenburg and Subramanian 2005a) pos- 
tulates that the deviations of the third-moment terms, 
TV'M*^^^^-' (k), from the contributions to these terms af- 
forded by the background turbulence, N'M^^^^'^ {\<i), are 
expressed through the similar deviations of the second 
moments, M(")(k) - M^"^(k): 

T{k) 

where T{k) is the characteristic relaxation time, which 
can be identified with the correlation time of the turbu- 
lent velocity field. The background turbulence is deter- 
mined by the budget equations and the general structure 
of the moments is obtained by symmetry reasoning. In 
the background turbulence, the mean magnetic field is 
zero. We applied the spectral r-approximation only for 
the nonhelical part hij of the tensor of magnetic fluctu- 
ations. We consider an intermediate nonlinearity which 
implies that the mean magnetic field is not enough strong 
in order to affect the correlation time of turbulent veloc- 
ity field. The theory for a very strong mean magnetic 
field can be corrected after taking into account a depen- 
dence of the correlation time of the turbulent velocity 
field on the mean magnetic field. 

We assume that the characteristic time of variation of 
the mean magnetic field B is substantially larger than the 
correlation time T{k) for all turbulence scales (which cor- 
responds to the mean-field approach). This allows us to 
get a stationary solution for the equations for the second 
moments fij , hij and gij . For the integration in k-space 
of these second moments we have to specify a model for 
the background turbulence (with B = 0) . We use a sim- 
ple model for the background homogeneous and isotropic 
turbulence. Using the derived equations for the second 
moments fij, hij and gij we calculate the mean electro- 
motive force £i = Eimn J 5rim(k)(ik. This procedure al- 
lows us to derive equations for the nonlinear coefficients 
defining the mean electromotive force in a homogeneous 
turbulence with a mean velocity shear (see for details, 
Rogachevskn and Kleeorin 2004, 2006). 

For simplicity in this study we do not take into account 
a quenching of the turbulent magnetic diffusion. This 
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FIG. 5: The dimensionless nonlinear coefficient a^^ (B) defin- 
ing tlie sfiear-current efTect for different values of the param- 
eter e: e = (solid); e = 1 (dashed). 

facet is discussed in details by Rogachevskii and Kleeorin 
(2004, 2006). We consider the nonlinear dynamo problem 
with the algebraic nonlinearity (B) which determines 
the nonhnear shear-current effect. The mean magnetic 
field is determined by the following nonlinear equations 

8 A 

— = Da,{B)B'^+A\ (13) 

^ = + K ' (14) 

(see Rogachevskii and Kleeorin 2004, Rogachevskii et al. 
2006), where B = |B|. The function a^{B) defining 
the nonlinear shear-current effect (which is normalized by 
(Tg ) , is shown in figure 5 for different values of the param- 
eter e. The asymptotic formulas for the nonlinear func- 
tion (T„ (B) are given by cr„ (B) — 1 for a weak mean mag- 
netic field {B < 5eq/4) andCT„(B) = -ll(H-e)/4(l-H9e) 
for B ^ Beq/4: (see figure 5). This implies that the non- 
linear function ct^^ (B) defining the shear-current effect 
changes its sign at some value of the mean magnetic field 
B = B^. Here B^ = 1.2Boq for e = and = 1.4Bcq 
for e = 1. However, there is no quenching of the non- 
linear shear-current effect contrary to the quenching of 
the nonlinear alpha effect, the nonlinear turbulent mag- 
netic diffusion, etc. The background magnetic fluctua- 
tions (caused by the small-scale dynamo and described 
by the parameter e), affect the nonlinear function cr„ (B). 

Numerical solutions of equations and H14|) for the 
nonlinear problem with the algebraic nonlinearity cr„ (B) 
are plotted in figures 6-8. In particular, these figures 
show the nonlinear evolution of the mean magnetic field 
B{t, z — Q) due to the shear-current effect for different 
values of the dynamo number D and the parameter e. 
The magnitude of the saturated mean magnetic field is 
in several times larger than the equipartition field de- 
pending on the dynamo number. Inspection of figures 7- 
8 shows that there is a range in the dynamo number 



FIG. 6: The nonlinear evolution of the mean magnetic field 
B{t, z = 0) due to the shear-current effect with the algebraic 
nonlinearity for e = (thin curve) and e = 1 (thick curve) 
and different near-threshold values of the dynamo number: 
D = 1.45 Dcr (solid); D = 1.3 Dcr (dashed); D = 1.15 Dcr 
(dashed-dotted). 
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FIG. 7: The nonlinear evolution of the mean magnetic field 
B{t, z = 0) due to the shear-current effect with the algebraic 
nonlinearity for e = (thin curve) and e = 1 (thick curve) 
and different values of the dynamo number: D — 10 (solid); 
7? = 30 (dashed); D = 50 (dashed-dotted). 

D — 22.8 — 59 where the nonlinear oscillations of mean 
magnetic field are observed at e = 0. 



4.2. The algebraic and dynamic nonlinearities 

In this study we consider nonhelical and nonrotating 
homogeneous turbulence. This implies that the kinetic 
helicity and the hydrodynamic a effect are zero. How- 
ever, the magnetic a effect caused by the small-scale 
magnetic helicity is not zero even in nonhelical turbu- 
lence. In particular, the magnetic helicity conservation 
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FIG. 8: The nonlinear evolution of the mean magnetic field 
B{t, 2 = 0) due to the shear-current effect with the algebraic 
nonlinearity for e = and different values of the dynamo 
number: (a). D ^ 47 (solid), D = 52 (dashed), D = 59 
(dashed-dotted); (b). D = 31 (solid), D = 35 (dashed); (c). 
D = 28 (solid); D = 22.77 (dashed); D = 22.765 (dashed- 
dotted). 



implies the growth of a magnetic alpha effect indepen- 
dent of whether kinetic helicity is driven into the system. 
In subsection 4.1 we have concentrated on the nonlinear 
shear-current effect (the algebraic nonlinearity) and have 
not discussed the effect of small-scale magnetic helicity 
(the dynamic nonlinearity) on the nonlinear saturation 
of the mean magnetic field. In this subsection we study 
joint action of the algebraic and dynamic nonlinearities. 

The small-scale magnetic helicity causes the magnetic 
a effect which is given by am = (i?) Xc(B), where 
$„(i?) is the quenching function of the magnetic a ef- 
fect given below by Eq. (|19|l . The function Xc(B) = 
(To/127r/9)(b-(Vxb)) is related to the small-scale cur- 
rent helicity (b-(Vxb)). For a weakly inhomogeneous 
turbulence (Iq <C L), the function Xc is proportional to 
the small-scale magnetic helicity, Xc = Xml {^^Trrj^p) (see 
Kleeorin and Rogachevskii 1999), where Xm = (s'b) is 
the small-scale magnetic helicity and a is the vector po- 
tential of small-scale magnetic field. The function Xc(B) 
entering the magnetic a effect is determined by the fol- 
lowing dynamical dimensionless equation (which is de- 
rived using arguments based on the magnetic helicity 
conservation law): 



dXc 
dt 



(15) 



(see, e.g., Kleeorin and Ruzmaikin 1982, Gruzinov and 
Diamond 1994, 1996, Kleeorin et al. 1995, 2000, 2002, 
2003a, 2003b, Kleeorin and Rogachevskii 1999, Blackman 
and Field 2000, Blackman and Brandenburg 2002, Bran- 
denburg and Subramanian 2005a, Zhang et al. 2006), 
where (l/3)(^o/-^)^Rin is the characteristic relax- 

ation time of the small-scale magnetic helicity, Rm is the 



magnetic Reynolds number and F is related to the flux 
of the small-scale magnetic helicity. 

The simplest form of the magnetic helicity flux is the 
turbulent diffusive flux of the magnetic helicity, F = 
— KyVxc (sec Kleeorin et al. 2002, 2003b), where the 
turbulent diffusivity coefficient measured in units 

of rj^ and the function Xc is measured in units of rjj,/ L. 
In real systems, the flux of small-scale magnetic helicity 
can be accompanied by some flux of large-scale magnetic 
helicity. However, the flux of large-scale magnetic he- 
licity does not explicitly enters in the dynamical equa- 
tion H15I) for the evolution of the small-scale magnetic 
helicity. This flux mostly affects the large-scale magnetic 
helicity. It can also introduce an additional anisotropy of 
turbulence, which can affect the dynamics of the mean 
magnetic field. 

Equation H15() determines the dynamics of the small- 
scale magnetic helicity, i.e. its production, dissipation 
and transport. For very large magnetic Reynolds num- 
bers (which are typical for many astrophysical situa- 
tions), the relaxation term Xcl^x very small, and it 
is very often dropped in Eq. H15|l in spite of the fact that 
the small yet finite magnetic diffusion is required for the 
reconnection of magnetic field lines. In particular, the 
magnetic Reynolds number, Rm does not enter into the 
steady state solution of Eq. (|15|) in the limit of very large 
Rm due to the effect of the magnetic helicity flux (see 
Kleeorin et al. 2003b). 

The account for the dynamics of the small-scale mag- 
netic helicity results in that the mean magnetic fleld is 
determined by the following dimensionless equations 



dA(t, z) 
dt~ 

dBy{t,z) 
dt 

dXc{t,z) 
dt 



^Da, (B) B' + By (B) Xc{t, z) + A" , 



-A' + B'^ , 



(16) 
(17) 



-n,x'^ + ^ = C[A' - By [D a, (B) B'^ 



^By^,{B)xc{t,z)- 



A"' 



(18) 



where am = ^ (B) Xc(i, is the magnetic a ef- 
fect and the parameter C — (2i/^o)^- In equations l|17|) 
and (UHJ we have neglected small terms ~ 0(5*^ ^) ^ 1. 
The algebraic function {B) in these equations is given 

by 



3 

8^2 



1 



arctan(\/8B) 



VSB 



(19) 



(see, e.g., Rogachevskii and Kleeorin 2000), where 
$„(B) = 1 - (24/5)^2 for B < l/VS and $„(B) = 
3/{8B'^) for B > l/%/8. Here the mean magnetic field 
B is measured in units of the equipartition field i?cq 
determined by the turbulent kinetic energy. In equa- 
tions (|16|I - H18|I there are four parameters: the dynamo 
number D — 'iS'^a^/C, the turbulent diffusivity of the 
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FIG. 9: The nonlinear evolution of the mean magnetic field 
B{t, 2 = 0) due to the shear-current effect with the algebraic 
and dynamic nonlinearities for e = 0, D = 2 Da, C — 100 
and difi'erent values of the parameter k^: Hj. — 0.5 (solid); 
— 0.3 (dashed); =0.1 (dashed-dotted). 



FIG. 11: The nonlinear evolution of the mean magnetic field 
B{t, 2 = 0) due to the shear-current effect with the algebraic 
and dynamic nonlinearities for e = 0, k^, = 10, D = 2 Da 
and different values of the parameter C: C = 100 (solid); 
C = 300 (dashed); C = 900 (dashed-dotted); C = 2700 
(dotted). 
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FIG. 10: The nonlinear evolution of the mean magnetic field 
B{t, z = 0) due to the shear-current effect with the algebraic 
and dynamic nonlinearities for e = 0, D — 2 Da, C = 100 
and difi^erent values of the parameter k^: = 10 (solid); 
hij. — 1 (dashed); =0.1 (dashed-dotted). 



small-scale magnetic helicity (measured in units of 
77y), the parameter C = (2L/lo)'^ and the relaxation time 
of the small-scale magnetic helicity — (4/3)Rm/C. 
Consider the simple boundary conditions for a layer of 



the thickness 2L in the z direction: By{t,\z\ — 1) 



0, 



A'{t,\z\ = 1) = and Xc(^, kl = 1) = 0, where z is 
measured in units of L. The initial conditions for the 
symmetric mode are chosen in the form: By{t = 0, z) = 
Bo cos(7rz/2), A{t = 0, z) = and Xcit = 0, z) = 0. 

Numerical solutions of equations H16|l - H18() for the non- 
linear problem with the algebraic and dynamic nonlinear- 
ities are plotted in figures 9-13. In particular, the non- 
linear evolution of the mean magnetic field B{t, z — 0) 



for different values of the parameters k^, , C, the dynamo 
numbers D and very large magnetic Reynolds numbers 
Km is shown in figures 9-13. Inspection of figures 9-10 
shows that the saturated level of the mean magnetic field 
depends strongly on the value of the turbulent diffusivity 
of the magnetic helicity k^,. The saturated level of the 
mean magnetic field changes from very small value for 

— 0.1 to the super-equipartition field for = 10. 
This is an indication of very important role of the trans- 
port of the magnetic helicity for the saturated level of 
the mean magnetic field. Indeed, during the generation 
of the mean magnetic field, the product 5-B is positive, 
and this produces negative contribution to the small-scale 
magnetic helicity (and negative magnetic a effect, see 
Eq. H15|) ). Therefore, this reduces the rate of genera- 
tion of large-scale magnetic field because the first and 
the second terms in the right hand side of Eq. H16() have 
opposite signs. The first term in Eq. H16|l describes the 
shear-current effect, while the second term in Eq. Hlti|) 
determines the magnetic a effect. If the magnetic helic- 
ity does not effectively transported out from the genera- 
tion region, the mean magnetic field is saturated even at 
small value of the magnetic field. Increase of the mag- 
netic helicity flux (by increasing of turbulent diffusivity 

of magnetic helicity) results in increase of the satu- 
rated level of the mean magnetic field above the equipar- 
tition field (see figures 9-10). Note that increase of the 
parameter C decreases the saturated level of the mean 
magnetic field (see figure 11). Actually the ratio k^/C 
determines the saturated level of the mean magnetic field 
in a steady-state. 

For the cases shown in figures 9-13, we drop the small 
relaxation term Xc/t^ in Eq. (|15|l due to very large mag- 
netic Reynolds numbers. Now we consider moderate 
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Table 1 

The saturated mean magnetic field 
versus the magnetic Reynolds number 



B/B 
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FIG. 12: The nonlinear evolution of the mean magnetic field 
B{t, 2 = 0) due to the shear-current effect with the algebraic 
and dynamic nonlinearities for e = 0, = 1, C = 100 and 
different values of the dynamo number; (a). D = 2 Da (thick 
solid); D = 1.7Dcr (dashed); D = 1.5Dcr (dashed-dotted); 
D = 1.43 Dcr (dotted); D = 1.1 (thin solid). 



magnetic Reynolds numbers, when the relaxation term 
Xc/tx in Eq. IjlSfl is not small but the flux of magnetic 
helicity is weak (e.g., Kj, = 0.1 — 0.3). In this case the 
small-scale magnetic helicity does not effectively trans- 
ported out from the generation region by the helicity 
flux. In Table 1 we demonstrate the effect of the moder- 
ate magnetic Reynolds numbers on the saturated level of 
the mean magnetic field. The decrease of the magnetic 
Reynolds numbers (and the relaxation time t^) increases 
the saturated level of the mean magnetic field (see Ta- 
ble 1), because the relaxation term Xc/tx in Eq. 115(1 de- 
creases the small-scale magnetic helicity in the generation 
region. On the other hand, for larger flux of small-scale 
magnetic helicity (k,, > 0.5), the effect of the magnetic 
Reynolds numbers on the saturated level of the mean 
magnetic field is very small. Note that the moderate 
magnetic Reynolds numbers Rm = 10 — 50 are irrelevant 
for astrophysical applications, although they are of an 
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FIG. 13: The nonlinear evolution of the mean magnetic field 
B(t, 2 = 0) due to the shear-current effect with the algebraic 
and dynamic nonlinearities for e = 0, = 1, C = 100 
and different values of the dynamo number: (a). D = 1 D^-c 
(solid); D = 5Dcr (dashed); D = 3-Dcr (dashed-dotted); 
D ^2D^^ (dotted). 



interest for the direct numerical simulations. 

Figures 12-13 show the nonlinear evolution of the mean 
magnetic field B{t, z = 0) for different values of the 
dynamo numbers D. The saturated level of the mean 
magnetic field increases with the increase of the dynamo 
numbers D within the range Z?cr < D < 2Dcr, and it 
decreases with the increase of the dynamo number for 
D > 2Dcr- This is a new feature in the nonlinear mean- 
field dynamo. For example, in the aQ dynamo the sat- 
urated level of the mean magnetic field usually increases 
with the increase the dynamo numbers. 

Generation of the large-scale magnetic field in a non- 
helical turbulence with an imposed mean velocity shear 
has been recently investigated by Brandenburg (2005) 
and Brandenburg et al. (2005) using direct numerical 
simulations. The results of these numerical simulations 
are in a good agreement with the numerical solutions of 
the nonlinear dynamo equations Hl()|l - H18|) discussed in 
section 4. 

Now let us compare the results of the numerical solu- 
tions of the nonlinear dynamo equations Hlt)l) - (|18|l with 
the numerical study by Brandenburg and Subramanian 
(2005b) of the mean-field dynamo with large-scale shear. 
In our study we use the expression for the nonlinear elec- 
tromotive force determined by Rogachevskii and Kleeorin 
(2004, 2006), which includes the nonlinear shear-current 
effect. On the other hand, Brandenburg and Subrama- 
nian (2005b) use very simplified form of the mean elec- 
tromotive force, neglecting e.g., the k effect related to the 
symmetric parts of the gradient tensor of the mean mag- 
netic field. This effect contributes to the shear-current 
effect (see Rogachevskii and Kleeorin 2003, 2004). Bran- 
denburg and Subramanian (2005b) have not taken into 
account the properties of the nonlinear shear-current ef- 
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feet found by Rogaehevskii and Kleeorin (2004, 2006). In 
partieular, there is no quenching of the nonlinear shear- 
current effect contrary to the quenching of the nonUn- 
ear alpha effect, the nonlinear turbulent magnetic diffu- 
sion, etc. During the nonlinear growth of the mean mag- 
netic field, the shear-current effect only changes its sign 
at some value of the mean magnetic field which affects 
the level of the saturated mean magnetic field (see Ro- 
gaehevskii and Kleeorin 2004, 2006). In our study we ne- 
glect small terms ~ 0(5""^) < 1 in equations ((TCl) - ((TH|l . 
i.e., we do not consider effect because the parameter 
S'* should be very large (see section 3). In addition we do 
not consider 5^ S effect because we neglected the small 
terms ^ 0[(Zo/-L)^] in equation ifTTji . Here the parameter 
S determines the WxJ term in the mean electromotive 
force. In our numerical solutions of the nonlinear mean- 
field dynamo equations we use the simplest form of the 
magnetic helicity flux (i.e., we use the turbulent diffusive 
flux of the magnetic helicity, F = — Vxcj where Krp IS 
not small), while Brandenburg and Subramanian (2005b) 
use the current helicity flux of Vishniac and Cho (2001) 
and a very small turbulent diffusive flux of the magnetic 
hehcity (with = 10~^). 

The parameter range in the study by Brandenburg and 
Subramanian (2005b) is different from that used in our 
study, and the maximum saturated level of the mean 
magnetic field obtained in the study by Brandenburg and 
Subramanian (2005b) is B = 0.6Boq at Rm = 10^ which 
strongly decreases with increase the magnetic Reynolds 
number Rm (see Table 4 of Brandenburg and Subra- 
manian 2005b). On the other hand, in our numeri- 
cal solutions of the nonlinear mean-field dynamo equa- 
tions which take into account the shear-current effect, 
the saturated level of the mean magnetic field reaches the 
super-equipartition field. These are the reasons why our 
numerical results discussed in this section are different 
from those obtained by Brandenburg and Subramanian 
(2005b). Note however, that increase of the saturated 
level of the mean magnetic field with the decrease of the 
magnetic Reynolds numbers in the case of very weak flux 
of small-scale magnetic helicity (found by Brandenburg 
and Subramanian 2005b), is confirmed by our numerical 
study (see Table 1). 



5. Discussion 

In this study we show that in a sheared nonhelical 
homogeneous turbulence the large-scale magnetic field 
can grow due to the shear-current effect from a very 
small seeding magnetic field. The shear-current dynamo 
strongly depends on the spatial scaling of the correla- 
tion time T{k) of the background turbulence. In partic- 
ular, for Kolmogorov scaling, T(fc) oc fc^^/^, the dynamo 
instability due to the shear-current effect occurs, while 
when r(fc) oc (for small hydrodynamic and magnetic 
Reynolds numbers) there is no the dynamo action in a 
sheared nonhelical turbulence. The dynamo instability 



is saturated by the nonlinear effects, and the dynami- 
cal nonlinearity due to the evolution of small-scale mag- 
netic helicity, plays a crucial role in nonlinear saturation 
of the large-scale magnetic field. The magnetic helicity 
flux strongly affects the saturated level of the mean mag- 
netic field in the nonlinear stage of the dynamo action. 
In particular, our numerical solutions of the nonlinear 
mean-field dynamo equations which take into account 
the shear-current effect, show that if the magnetic he- 
licity flux is not small, the saturated level of the mean 
magnetic field is of the order of the equipartition field 
determined by the turbulent kinetic energy. 

The shear-current dynamo acts also in inhomogeneous 
turbulence. However, in inhomogeneous turbulence with 
a large-scale velocity shear the kinetic helicity and hy- 
drodynamic a effect are not zero (see Rogaehevskii and 
Kleeorin 2003; 2006; Radler and Stepanov 2006). In this 
case the shear-current dynamo acts together with the a- 
shear dynamo (which is similar to the afl dynamo). The 
joint action of the shear-current and the a-shear dynamo 
has been recently discussed by Rogaehevskii and Kleeorin 
(2006); Pipin (2006). 

Turbulence with a large-scale velocity shear is a uni- 
versal feature in astrophysics, and the obtained results 
can be important for explanation of the large-scale mag- 
netic fields generated in astrophysical sheared turbu- 
lence. Rogaehevskii et al. (2006) have suggested that 
the shear-current effect might be considered as an origin 
for the large-scale magnetic fields in colliding protogalac- 
tic clouds and in merging protostellar clouds. 

Note that the problem of the generation of the mean 
magnetic field in a turbulence with large-scale velocity 
shear is similar to that for generation of mean vortic- 
ity in a sheared turbulence. The instability of the per- 
turbations of the mean vorticity in a turbulence with a 
large-scale linear velocity shear was studied by Elperin 
et al. (2003). This instability is caused by a combined 
effect of the large-scale shear motions (skew-induced de- 
flection of equilibrium mean vorticity due to the shear) 
and Reynolds-stress-induced generation of perturbations 
of mean vorticity. In Appendix A we compare the prob- 
lem of generation of the mean magnetic field in a turbu- 
lence with a large-scale velocity shear with that of the 
generation of mean vorticity in a sheared turbulence. 
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APPENDIX A: THRESHOLD FOR GENERATION 
OF MEAN VORTICITY IN A SHEARED 
TURBULENCE 

The problem of generation of the mean magnetic field 
in a turbulence with large-scale velocity shear is simi- 
lar to that of generation of mean vorticity in a sheared 
turbulence. Indeed, let us discuss the generation of the 
mean vorticity in a turbulence with a large-scale linear 
velocity shear, U = (0, Sx, 0) and W = (0, 0, S). Pertur- 
bations of the mean vorticity W = {Wx{z),Wy{z),0) are 
determined by the following equations: 



dt 



dt 



y = -a,,SllW': + v,W'' 



(A.l) 
(A.2) 



(see for details Elperin et al. 2003), where v.^ is the tur- 
bulent viscosity, W^' = d^Wi/dz^ and the parameter a„ 
is given by Eq. (|A.4|I below. The first term, SWy, in the 
right hand side of Eq. I|A.1(I determines a skew-induced 
generation of perturbations of the mean vorticity Wx by 
deflection of the equilibrium mean vorticity W, where 
U are the perturbations of the mean velocity. In partic- 
ular, the mean vorticity Wx^x is generated from WyBy 
by equilibrium shear motions with the mean vorticity 
W ^ Se^, i.e., SWyBx cx (W-V)?7:,e:, cx WyBy x W. 
Here ex, By and are the unit vectors along x, y and 
z axis, respectively. On the other hand, the flrst term, 
-o-„, SIq Wx, in the right hand side of Eq. I|A.2|I deter- 
mines a Reynolds-stress-induced generation of perturba- 
tions of the mean vorticity Wy by turbulent Reynolds 
stresses. This implies that the mean vorticity WyBy 
is generated by an effective anisotropic viscous term 
oc -llA{WxBx-V)U{x)ey cx -l^SWi!By. This mech- 
anism of the generation of perturbations of the mean 
vorticity WyBy can be interpreted as a stretching of the 
perturbations of the mean vorticity WxBx by the equilib- 
rium shear motions U = S xBy during the turnover time 
of turbulent eddies (see Elperin et al. 2003). 

Note that equations (|A.l|) - (|A.2p for the perturbations 
of the mean vorticity are very similar to equations (Q- 
^ for the perturbations of the mean magnetic fleld in a 
sheared turbulence. The growth rate 7 of the instability 



13 



of the perturbations of the mean vorticity is given by 



(A.3) 



The form of the growth rate I) A. 3(1 of the perturbations 
of the mean vorticity is also very similar to the growth 
rate (O of the mean magnetic field due to the shear- 
current effect. On the other hand, the magnetic dynamo 
instability is different from the instability of the pertur- 
bations of the mean vorticity although they are governed 
by similar equations. The mean vorticity W = VxU 
is directly determined by the velocity field U, while the 
magnetic field depends on the velocity field through the 
induction equation and Navier-Stokes equation. 

In the present study we derived a more general form 
of the parameter cr^ for an arbitrary scaling of the cor- 
relation time T{k) of the turbulent velocity field. The 
parameter cr„, is derived using Eq. (21) of the paper by 
Elperin et al. (2003). It is given by 



4?o 
45 



^0 



43- 



63 



(A.4) 



where / and /q are determined by Eq. |0J. The instability Eq. lO] . 



depends on the correlation time T(fc). In particular, when 
T(fc) oc fc"'^, the ratio ///q ~ — /i, and the criterion of the 
instability reads 2 /i'^ — 43 /i + 63 > 0, i.e., the instability 
is excited when < ^ < 1.58 and /i > 19.9. Note that 
the condition ^ > 19.9 is not realistic. 

When r(/c) cx fc^"', we recover the result obtained by 
Elperin et al. (2003), i.e., = 4(2g2 _ 47g 108)/315. 
In particular, for the Kolmogorov scaling, T(fc) cx fc^^/"^ 
(i.e., for q — 5/3), we arrive at « 0.45. 

For small hydrodynamic Reynolds numbers, the scal- 
ing T{k) - l/ivk"^), the ratio I/Io = -2 (i.e., ^ = 2), 
and the parameter = —4/9 < 0. This implies that 
the instability of the perturbations of the mean vorticity 
does not occur in agreement with the recent results by 
Riidiger and Kichatinov (2006). They have not found the 
instability of the perturbations of the mean vorticity in 
a random flow with large-scale velocity shear using the 
second order correlation approximation (SOCA). Note 
that this approximation is valid only for small hydrody- 
namic Reynolds numbers [see discussion in section 3 after 



